Observation of the antimatter lieiium-4 nucieus ( submitted to Na- 
ture, under media embargo) 

High-energy nuclear collisions create an energy density similar to that of the universe mi- 
croseconds after the Big Bang^, and in both cases, matter and antimatter are formed with 
comparable abundance. However, the relatively short-lived expansion in nuclear collisions 
allows antimatter to decouple quickly from matter, and avoid annihilation. Thus, a high- 
energy accelerator of heavy nuclei is an efficient means of producing and studying antimatter. 
The antimatter helium-4 nucleus C^He), also known as the anti-a (a), consists of two antipro- 
tons and two antineutrons (baryon number B = —4). It has not been observed previously, 
although the a particle was identified a century ago by Rutherford and is present in cosmic 
radiation at the 10% leveF. Antimatter nuclei with B < —1 have been observed only as 
rare products of interactions at particle accelerators, where the rate of antinucleus produc- 
tion in high-energy collisions decreases by about 1000 with each additional antinucleoiP^. 
We present the observation of the antimatter helium-4 nucleus, the heaviest observed antinu- 
cleus. In total 18 ^He counts were detected at the STAR experiment at RHld^in 10^ recorded 
gold on gold (Au+Au) collisions at center-of-mass energies of 200 GeV and 62 GeV per 
nucleon-nucleon pair. The yield is consistent with expectations from thermodynamic' and 
coalescent nucleosynthesis^ models, which has implications beyond nuclear physics. 

In 1928, DiracP predicted the existence of negative energy states of electrons based on the 
application of symmetry principles to quantum mechanics, but only recognized these states to be 
antimatter after Anderson'^ discovery of the positron (the antielectron) in cosmic radiation four 
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years later. The predicted antiprotons'^ and antineutronJ^^ were observed in 1955, followed by 
antideuterons (d), antitritons ('^H), and antihelium-3 ("^He) during the following two decade^^^^. 
Recent accelerator and detector advances led to the first production of antihydrogen'^ atoms in 
1995 and the discovery of strange antimatter, the antihypertriton (|H), in 2010 at the Relativistic 
Heavy Ion Collider (RHIC) at Brookhaven National Laboratory (BNL^. 

Collisions of relativistic heavy nuclei create suitable conditions for producing antinuclei, be- 
cause large amounts of energy are deposited into a more extended volume^ than that achieved in 
elementary particle collisions. These nuclear interactions briefly (~ 10"^^ seconds) produce hot 
and dense matter containing roughly equal numbers of quarks and antiquark^, often interpreted 
as quark gluon plasma (QGPj^. In contrast to the Big Bang, nuclear collisions produce negli- 
gible gravitational attraction and allow the plasma to expand rapidly. The hot and dense matter 
cools down and transitions into a hadron gas, producing nucleons and their antiparticles. The pro- 
duction of light antinuclei can be modeled successfully by macroscopic thermodynamic^, which 
assumes energy equipartition, or by a microscopic coalescence procesd^^, which assumes uncor- 
related probabilities for antinucleons close in position and momentum to become bound. The high 
temperature and high antibaryon density of relativistic heavy ion collisions provide a favorable 
environment for both production mechanisms. 

The central detector used in our measurements of antimatter, the Time Projection Chamber 
(TPCj^of the STAR experiment (Solenoidal Tracker At RHIC), is situated in a solenoidal mag- 
netic field and is used for three-dimensional imaging of the ionization trail left along the path of 
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Figure 1 : A three-dimensional rendering of the STAR TPC surrounded by the TOF barrel shown as 
the outermost cylinder. Tracks from an event which contains a ^He are shown, with the ^He track 
highlighted in bold red. 

charged particles (Fig. [T]). In addition to the momentum provided by the track curvature in the 
magnetic field, the detection of ''He particles relies on two key measurements: the mean energy 
loss per unit track length (dE/dx) in the TPC gas, which helps distinguish particles with different 
masses or charges, and the time of flight of particles arriving at the time of flight barrel (TOFj^^l 
surrounding the TPC. In general, time of flight provides particle identification in a higher momen- 
tum range than (dE/dx). The (dE/dx) resolution is 7.5% and the timing resolution for TOF is 95 
picoseconds. 

The trigger system at STAR selects collisions of interest for analysis. The minimum-bias 



3 



(MB) trigger selects all particle-producing collisions, regardless of the extent of overlap of the 
incident nuclei. A central trigger (CENT) preferentially selects head-on collisions, rejecting about 
90% of the events acquired using the MB trigger. The sample of 10^ Au-i-Au collisions used 
in this search is selected based on MB, CENT, and on various specialized triggers. Preferential 
selection of events containing tracks with charge Ze = ±2e (where e is the electron charge) 
was implemented using a High-Level Trigger (HLT) for data acquired in 2010. The HLT used 
computational resources at STAR to perform a real-time fast track reconstruction to tag events that 
had at least one track with a (dE/dx) value that is larger than a threshold set to three standard 
deviations below the theoretically expected valueP^ for ^He at the same momentum. The HLT 
successfully identified 70% of the events where a "^He track was present while selecting only 0.4% 
of the events for express analyses. 
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Figure 2: {dE/dx) versus p/\Z\ for negatively charged particles (left) and positively charged 
particles (right). The black curves show the expected values for each species. The lower edges 
of the colored bands correspond to the HLT's online calculation of 3cr below the (dE/dx) band 
centei^^ for •^He. 

Figure [2] shows (dE/dx) versus the magnitude of magnetic rigidity, p/\Z\, where p is mo- 
mentum. A distinct band of positive particles centered around the expected valu^for ^He particles 
indicates that the detector is well-calibrated. On the left side of Fig.|2} where p/\Z\ is less than 1.4 
GeV/c, four negative particles are particularly well separated from the '^He band and are located 
within the expected band for ^He. Above 1.75 GeV/c, (dE/dx) values of ^He and ^He merge 
and the TOF system is needed to separate these two species. 

The top two panels of Fig. [3] shows the (dE/dx) (in units of multiples of cXdE/dx, iT-a-aE/dJ 
versus calculated mass m = (p/ c)^ {t?c^ / L"^ — 1), where adE/dx is the rms width of the (dE/dx) 
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distribution for ^He or ^He, t and L are the time of flight and path length, respectively, and c is 
the speed of light. The first (second) panel shows negatively (positively) charged particles. In both 
panels, majority species are ^He and ^He. In the second panel, the ^He particles cluster around 
'^'^dE/dx ^ ^^'^ mass = 3.73 GeV/c^, the appropriate mass for ^He. A similar but smaller cluster 
of particles can be found in the first panel for ^He. The bottom panel shows the projection onto 
the mass axis of the top two panels for particles with n„^^^^^^ of -2 to 3. There is clear separation 
between ^He and ^He mass peaks. Eighteen counts for ^He are observed. Of those, sixteen are 
from collisions recorded in 2010. Two counts^^ identified by (dE/dx) alone from data recorded in 
2007 are not included in this figure, because the STAR TOP was not installed at that time. 
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Figure 3: The top two panels show the {dE jdx) in units of multiples of adE/dx, fT-a^iE/dx' °f '^^S" 
atively charged particles (first panel) and positively charged particles (second panel) as a func- 
tion of mass measured by the TOF system. The masses of '^He (^He) and "^He (^He) are indi- 
cated by the vertical lines at 2.81 GeV/c^ and 3.73 GeV/c^, respectively. The horizontal line 
marks the position of zero deviation from the expected value of (dE/dx) {n^^^^^^^ = 0) for 
^He (^He). The rectangular boxes highlight areas for ''He C'He) selections : —2 < n^^^^^^^ < 3 
and 3.35 GeV/c^ < mass < 4.04 GeV/c^ (corresponding to a ±3(t window in mass). The bottom 
panel shows a projection of entries in the upper two panels onto the mass axis for particles in the 
window of —2 < o dE/dx < 3. The combined measurements of energy loss and the time of flight 
allow a clean identification to be made in a sample of 0.5 x 10^^ tracks from 10^ Au-i-Au collisions. 



To evaluate the background in ^He due to ^He contamination, we simulate the ^He mass dis- 
tribution with momenta and path lengths, as well as the expected time of flight from '^He particles 
with timing resolution derived from the same data sample. The contamination from misidentifying 
^He as ^He is estimated by integrating over the region of the ^He selection. We estimate that the 
background contributes 1.4 (0.05) counts of the 15 (1) total counts from Au+Au collisions at 200 
(62) GeV recorded in 2010. 

The observed counts are used to calculate the antimatter yield with appropriate normaliza- 
tion (the differential invariant yield) in order to compare to the theoretical expectation. Detec- 
tor acceptance, efficiency, and antimatter annihilation with the detector material are taken into 
account when computing yields. Various uncertainties related to tracking in the TPC, match- 
ing in the TOF, and triggering in the HLT are cancelled when the yield ratios of '^He/'^He and 
^He/^He are calculated. The ratios are ^He/^He = (3.0 ± 1.3(stat)t[j|(sys)) x lO^^ 
^H^/^H^ = (3.2 ± 2.3(stat)i[ji(sys)) x IQ-^ for central Au-i-Au collisions at 200 GeV. The 
ratios were obtained in a window 40° < 9 < 140°, where polar angle, 9, is the angle between the 
particle's momentum vector and the beam axis (these 9 limits correspond to limits of —1 to 1 in a 
related quantity, pseudorapidity), and in ap-r per baryon window centered atpT/\B\ = 0.875 GeV/c 
with a width of 0.25 GeV/c, where pr is the projection of the momentum vector on the plane that 
is transverse to the beam axis. Ratios calculated by a Blastwave modeP^for the pt/\B\ window 
mentioned above and for the whole range of pt/\B\ differ by only 1%. The differential yields 
{(PN/ (pTdpTdy)) for ^He (^He) are obtained by multiplying the ratio of '^He/'^He (^He/^He) with 
the ^He (^He) yieldd^. The systematic uncertainties consist of background (-6% for both ratios), 
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feed-down from (anti-)hypertritons (18% for both ^He and ^He), knock-outs from beam-material 
interactions (-5% for the ratio ^He/^He only) and absorption (4% for the ratio ''He/'^He only). Fig- 
ure [4] shows the exponentiaP invariant yields versus baryon number in 200 GeV central Au-i-Au 
collisions. Empirically, the production rate reduces by a factor of l.Gloie ^ (l-l^oi ^ 10'^) 
for each additional antinucleon (nucleon) added to the antinucleus (nucleus). This general trend is 
expected from coalescent nucleosynthesis models^, originally developed to describe production of 
antideuteronP^, and as well as from thermodynamic modelP. 
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Figure 4: Differential invariant yields as a function of baryon number B, evaluated at pt/\B\ = 
0.875 GeV/c, in central 200 GeV Au-i-Au collisions. Yields for (anti)tritons (^H and ^H) lie close 
to the positions for '^He and ^He, but are not included here because of poorer identification of 
(anti)tritons. The lines represent fits with the exponential formula oc e"'''^' for positive and nega- 
tive particles separately, where r is the production reduction factor. Analysis details of yields other 
than ''He C'He) have been presented elsewherd^HH. Errors are statistical only. Systematic errors are 
smaller than the symbol size, and are not plotted. 
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In a microscopic picture, a light nucleus emerging from a relativistic heavy-ion collision 
is produced during the last stage of the collision process. The quantum wave functions of the 
constituent nucleons, if close enough in momentum and coordinate space, will overlap to produce 
the nucleus. The production rate for a nucleus with baryon number B is proportional to the nucleon 
density in momentum and coordinate space, raised to the power of \B\, and therefore exhibits 
exponential behaviour as a function of B. Alternatively, in a thermodynamic model, a nucleus is 
regarded as an object with energy _E ~ \B\mj\f, where m^v is the nucleon mass, and the production 
rate is determined by the Boltzmann factor exp {—E/T), where T is the temperature^Il. This 
model also produces an exponential yield. A more rigorous calculation!^ can provide a good fit to 
the available particle yields and predicts the ratios integrated over px to be "^He/^He = 3.1 x 10^^ 
and ^He/^He = 2.4 x lO^'^, consistent with our measurements. The considerations outlined above 
offer a good estimate for the production rate of even heavier antinuclei. For example, the yield of 
the stable antimatter nucleus next in line (B = —6) is predicted to be down by a factor of 2.6 x 10^ 
compared to ^He and is beyond the reach of current accelerator technology. 

A potentially more copious production mechanism for heavier antimatter is by the direct 
excitation of complex nuclear structures from the vacuum^^. A deviation from the usual rate re- 
duction with increasing mass would be an indication of a radically new production mechanisnP. 
On the other hand, going beyond nuclear physics, the sensitivity of current and planned space- 
based charged particle detectors is below what would be needed to observe antihelium produced 
by nuclear interactions in the cosmos, and consequently, any observation of antihelium or even 
heavier antinuclei in space would indicate the existence of a large amount of antimatter elsewhere 
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in the universe. In particular, finding ^He in the cosmos is one of the major motivations for space 
detectors such as the Alpha Magnetic Spectrometei^. We have shown that ''He exists, and have 
measured its rate of production in nuclear interactions, providing a point of reference for possi- 
ble future observations in cosmic radiation. Barring one of those dramatic discoveries mentioned 
above or new breakthrough in the accelerator technology, it is likely that the ''He will remain the 
heaviest stable antimatter nucleus observed for the foreseeable future. 
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